Flow reversals in thermally driven turbulence 
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We analyze the reversals of the large scale flow in Rayleigh-Benard convection both through par- 
ticle image velocimetry flow visualization and direct numerical simulations (DNS) of the underlying 
Boussinesq equations in a (quasi) two-dimensional, rectangular geometry of aspect ratio 1. For 
medium Prandtl number there is a diagonal large scale convection roll and two smaller secondary 
rolls in the two remaining corners diagonally opposing each other. These corner flow rolls play a 
crucial role for the large scale wind reversal: They grow in kinetic energy and thus also in size thanks 
to plume detachments from the boundary layers up to the time that they take over the main, large 
scale diagonal flow, thus leading to reversal. Based on this mechanism we identify a typical time 
scale for the reversals. We map out the Rayleigh number vs Prandtl number phase space and find 
that the occurrence of reversals very sensitively depends on these parameters. 



Spontaneous flow reversals occur in various buoyancy 
driven fluid dynamical systems, including large scale 
flows in the ocean, the atmosphere, or the inner core 
of stars or the earth, where such reversals are associated 
with the reversal of the magnetic field. The paradig- 
matic example for buoyancy driven flow is the Rayleigh- 
Benard system, i.e., a fluid- filled cell heated from below 
and cooled from above, see e.g. the recent reviews [T, "2^. 
In this system flow reversals have been detected and sta- 
tistically analyzed, mainly through measurements of the 
temperature at one or several points [T, ^ in the flow 
or at the walls and more recently through flow visualiza- 
tion with particle image velocimetry (PIV) [6l[7]. Various 
models have been developed to account for the reversals, 
either of stochastic nature [HI [9] or based on simplify- 
ing (nonlinear) dynamical equations [TOl [TT] , which show 
chaotic deterministic behavior. Most of the experimental 
studies have so far been done in a cylindrical cell, where 
the complicated three-dimensional dynamics of the con- 
vection role (see e.g. [12j [13] and section VIII of [Ij) 
complicates the identification of the reversal process. 

In the present paper, we restrict ourselves to the study 
of flow reversals in (quasi) two-dimensional (2D) rectan- 
gular geometry: experimentally to RB convection in a 
flat cell and numerically to DNS of the two-dimensional 
Oberbeck-Boussinesq equations, for which reversals have 
been observed already in [14] . This approach offers three 
advantages: (i) The flow reversal in quasi-2D is less com- 
plicated than in 3D (and therefore of course may be dif- 
ferent); (ii) the visualization of the full flow is possible; 
and (iii) a study of a considerable fraction of the Rayleigh 
number Ra vs Prandtl number Pr phase space becomes 
numerically feasible. 

The experiments were made in rectangular, quasi-2D 
cells [TF . To extend the range of Ra, two cells of nearly 



identical geometry are used. The larger (smaller) cell has 
a horizontal cross section of 24.8x7.5 (12.6x3.8) cm^, 
and the heights of the larger ( smaller) cell is H = 25 A 
cm {H = 12.6 cm), giving an aspect ratio F ^ 1 in the 
plane of the main flow (and an aspect ratio of about 0.3 
in the direction perpendicular to it). The fluid is wa- 
ter with a mean temperature of 28^ C, corresponding to 
Pr=5.7. For direct visualization of flow reversals, PIV 
is used for a few selected Ra. The PIV measurement 
in this system has also been described previously [T5] . 
To study the statistical properties of the reversals over 
long time periods, we measure the temperature contrast 
5T between two thermal probes embedded respectively 
in the left and right sides of either the top or the bottom 
plates. Reversals of the upward going hot plumes and 
downward going cold plumes correspond to the switch- 
ing between the right and left sides of the system, ST 
therefore is indicative of reversals. 

The numerical code is based on a fourth-order finite- 
difference discretization of the incompressible Oberbeck- 
Bousinesq equations and has been described in [16] . The 
grid resolution has been chosen according to the strict 
requirements as formulated in \17^. As in experiment 
also the numerical flow is wall-bounded, i.e., we use no- 
slip boundary conditions at all solid boundaries: u = 
at the top {z = H) and bottom {z = 0) plates as well as 
on the side walls x = and x = H. For the temperature 
at the side walls heat-insulating conditions are employed 
and T5 — = A is the temperature drop across the 
whole cell. Times are given in multiples of the large eddy 
turnover time t^;, defined by Ie 47r/ {\ujc{t)\)^ where 
ujc denotes the center vorticity. 

We start by showing qualitative features of the reversal 
process using examples from both numerical simulations 
and experiments. FigjT] shows snapshots of the tempera- 
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FIG. 1. (Color online) Top panel: Snapshots of the temperature (color) and velocity (arrows) field and time trace of angular 
momentum from numerical simulations (Ra — 10^ and Pr — 4.3). Bottom panel: Snapshots of the velocity field and time 
trace of angular momentum from experiment {Ra — 3.8 x 10^ and Pr — 5.7). (a), (b), and (c) show the instantaneous 
dimensionless temperature (T — Tt) / /S. distribution, (d) shows the temporal change of the dimensionless angular momentum 
L(t)/Lo, where Lq is the maximum of the absolute value of L. The positive and negative signs indicate the anti-clockwise and 
clockwise circulations, respectively, (e), (f), and (g) show the PIV-measured instantaneous velocity field and (h) the normalized 
instantaneous angular momentum L{t)/Lo. The color bar indicates the magnitude of the velocity (in unit of cm/sec). The 
snapshots in (a), (b), and (c) give numerical, and those in (e), (f), and (g) give experimental examples of the large scale 
circulation before, during and after a reversal process, as indicated in (d) and (h) respectively. Note that (b) and (f) show 
clearly the key role played by the growth of the corner rolls in the reversal process. Corresponding movies are offered in the 
supplementary material. 



ture and velocity fields from DNS and those of the veloc- 
ity field from experiment just before, during, and after 
the large convection roll reversal. Corresponding videos 
can be viewed from the supplementary materials. Vi- 
sually, the reversal process can be easily detected. To 
automatize this we measure the local angular velocity at 
the center of the cell; however, with this method some 
plumes passing through the center can lead to erroneous 
reversal counting. A better way is to rely on a global 
quantity, e. g. the global angular momentum (which 
has been successfully used for reversal characterization 
in 2D Navier-Stokes turbulence ^S1[18|). This is defined 
as L{t) = {-{z - H/2)u^{x,t) + (x - H/2)u,{x,t))v, 
where represents averaging over the full volume. 

The time dependence of L from simulation and experi- 
ment, as shown respectively in Figsjljd) and (h), indeed 
nicely reveals the reversal through a sign-change. 

From the movies corresponding to Fig. [l] the basic role 
of the corner flows in the reversal process can be ob- 
served: While the main roll is diagonally orientated in 
the cell, smaller counter-rotating rolls develop in diag- 
onally opposing corners. They are energetically fed by 
detaching plumes from the boundary layers trapped in 
the corner flows, leading to their growth. Once the two 
corner flows have reached a linear extension of ~ H/2 



[FigjTJb) and moment (b) in FigjTJd) and Fig. [TJf) and 
moment (f) in Figjljh)], they destroy the main convec- 
tion roll and establish another one circulating in opposite 
direction. 

The heights h{t) of the corner flows are measured by 
first fitting the temperature profile at the respective side- 
wall with splines, and then identifying the position of the 
steepest gradient of T{z): From movies and snapshots 
we judge that this is an excellent measure for the height 
h{t) of the corner- flow. Time series of /i(t), together with 
the (rescaled) center vorticity uJc{t) as quantitative mea- 
sure of the strength of the large scale convection roll, are 
shown in figure |2j It is seen that after a reversal (as indi- 
cated by a sign change in 0Jc{t)) the respective corner-flow 
grows roughly linearly in time, before it reaches the half- 
height h{t)/H ^1/2 and breaks down, leading to flow 
reversal. However, the growth of the corner-flow need 
not always lead to a reversal of the large scale convec- 
tion roll: There are cases in which the corner-flow looses 
energy due to some plume detachment from it, leading 
to full recovery of the large scale convection roll in its 
orginal direction (e.g., at t/tE ~ 1300 in figure |2|. Also 
in experiment we have observed such unsuccessful build- 
ups of the corner-flow. Below we will try to quantify the 
energy gains and losses of the corner- flows. 
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FIG. 2. (Color online) Time series of the center vorticity uJc{t) 
(rescaled by its maximum) (lower panel) and the heights h(t) 
of the lower left (blue) and right (red) corner flows, revealing 
their approximate linear growth. Not all growth processes 
need to lead to an immediate successful reversal, as seen for 
t/tE ^ 1300. Ra 10^ Pr 4.3. 



FIG. 3. (Color online) Log- log plot of the Ra-dependence 
of the mean time intervals (r) between reversals, normalized 
in terms of the large eddy turnover time Ie- Red symbols 
are from experiment and black ones from simulation. The 
error bar originates from the statistics of the reversals; for 
the numerical case it is smaller than the symbol size. 



The mean time interval (r) is shown in Fig. [s] Firstly, 
we clearly see that experiment and simulation are in very 
good agreement. The figure shows that (r) /tE at most 
weakly depends on Ra up to Ra ^ 2 x 10^, but for larger 
Ra rapidly increases with increasing Ra, i.e., reversals 
occur less and less frequently. The numbers mean that 
there are only very few reversals: of order one per hour 
in the Ra ~ 10^ range down to one within two days in 
the Ra ~ 10^ range. For Ra> 5 x 10^ no reversals could 
be detected any more in our numerical simulations, even 
for an averaging time of 2600 large eddy turnovers (see 
accompanying movie). In experiment two reversals could 
still be observed at Ra= 1.6 x 10^ (presumably due to the 
longer observation time in experiment which goes beyond 
10000 large eddy turnovers, corresponding to four days), 
but these also cease for larger Ra. 

These findings led us to map out a considerable frac- 
tion of the Ra-Pr parameter space. The result is shown 
in figure [4j where the black symbols are from simulation 
and red symbols are from experiment. One sees a rather 
complicated structure. But given the limited amount of 
data, experiment and simulation are in general agree- 
ment, especially considering the fact that the simulations 
are for the true 2D case whereas the experiments run in a 
quasi-2D cell. It should again be pointed out that the ex- 
perimental data point with the highest Ra (= 1.6 x 10^) 
that still shows a reversal has an extremely low reversal 
rate (0.5/day), which corresponds to waiting for ~ 5500 
large scale turn over time for a single reversal to occur. 

From figure [4] we conclude that not only for too large 
Ra (as compared to above case of figure 1 with Ra = 10^ 
and Pr = 4.3) the reversals do not occur any more, but 
also for too large or too small Pr. How to physically un- 
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FIG. 4. (Color online) Phase diagrams in the Ra-Pr plane. 
Red symbols are from experiment and the black ones from 
DNS. Circles correspond to detected reversals [23], crosses 
to no detected reversals, in spite of excessive simulation (or 
observation) time. The straight lines are guides to the eye; 
they have (from left to right) slopes 0.25 and 1.00. 



derstand this complicated behavior? The key towards an 
understanding lies, from our point of view, in the role of 
the corner flows, and is based on a detailed observation 
of many movies at various Ra and Pr (see accompany- 
ing material). As stated above, the corner flow rolls are 
fed by plumes detaching from the plates' boundary lay- 
ers. For too small Pr (i.e., too large thermal diffusivity) 
the thermal energy they carry is lost through thermal 
diffusion. On the other hand, for too large Pr (i.e., too 
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large kinematic viscosity) the thermal BL is nested in 
the kinematic BL and the thermal coupling of the cor- 
ner flow towards the thermal BL is hindered. In both 
cases the buildup of the corner flow and thus the rever- 
sals are suppressed. The situation is similar to the one in 
rotating RB, where Ekman vortices form, sucking heat 
out of the thermal BL and enhancing the heat flux: Also 
here there is an optimal Pr ~ 10 for which the Nusselt 
number is maximal, and for larger or smaller Pr the very 
same above mechanisms hinder efficient heat transport 

We now quantify this argument. The heat influx feed- 
ing the corner flow scales as Jin ~ n/S.H~^Nu. The out- 
flux of thermal and kinetic energy is either of diffusive or 
of convective origin. We model it as J out = Jtut ^ ^mlt^ - 
Flow reversal is prevented if J out > Jin- The convective 
outfiux, which is dominant for large Pr, is modeled by 
Jour - Jin ■ K/Xe - K./^H-^Nu^/^/Re. The diffusive 
outfiux is JquJ ^ nt/S.H~^Nu with some effective, tur- 
bulent thermal diffusivity Kt = I't/Pr ~ Pr~^U^/e ~ 
nPr'^Re^ / {NuRa)^ where we have assumed Pvt ~ Pr 
and employed the /c-e-model [20 for the turbulent vis- 
cosity Vf 

For dominant diffusive outfiux (thus low Pr), suppres- 
sion of reversals occurs at nt > n. The threshold is 
given by the scaling relation NuRa ~ Pr^Re^. Insert- 
ing Nu{Ra, Pr) and Re{Ra, Pr) either from experiment 
or from the unifying theory of refs. [21 , one obtains a 
relation between the critical Prandtl number Prcrit and 
the critical Rayleigh number Racrit at which reversals 
stop. Depending on whether regime IIu, or IVu of 
refs. [21] is dominant, we obtain Prcrit ^ J^^Zrit^ ^iih. 
7 = 3/5 or 2/3, respectively, which correctly reflects the 
trend in Fig. [4j 

For large Pr the convective outfiux will be dominant. 
Here the threshold condition is ~ nNu/ \/Re^ which 
with Nu{Ra^ Pr) and Re{Ra^ Pr) in regime lu of refs. 
[21 leads to an Ra-independent Prcrit^ beyond which no 
reversals are possible. The reality of figure [4] is clearly 
more complicated [22 but at least the general trends are 
consistent with this explanation. 

Finally, we note that we also performed experiments 
and simulations for F = 0.85. Even for this relatively 
small change in F the overall flow dynamics is very dif- 
ferent and much more complex as compared to the case 
of F = 1. Just as the important role the corner flows play 
for reversals, this finding demonstrates the strong effect 
of the cell geometry on the overall flow dynamics in the 
F = 0{oo) regime. In full 3D geometries, it may be less 



pronounced, but it certainly is present, too, see also ref. 
[23\ It remains remarkably that the rich structure in 
the (Ra, Pr, F) parameter space for reversals is hardly 
reflected in Nu and Re. 
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